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ABSTRACT 

We report optical Cousins R and I band monitoring observations of the high redshift (z 
= 4.67) QSO SDSS J153259.96-003944.1 that does not show detectable emission lines 
in its optical spectrum. We show this object varies with a maximum amplitude of ^0.4 
mag during a year and three months of monitoring. Combined with two other epochs 
of photometric data available in the literature, we show the object has gradually faded 
by ~0.9 mag during the period June 1998 — April 2001. A linear least squares fit to all 
available observations gives a slope of ~0.35 mag/yr which translates to ~1.9 mag/yr 
in the rest frame of the quasar. Such a variability is higher than that typically seen 
in QSOs but consistent with that of BL Lacs, suggesting that the optical continuum 
is Doppler boosted. Alternatively, within photometric errors, the observed lightcurve 
is also consistent with the object going through a microlensing event. Photoionization 
model calculations show the mass of the Broad Line Region to be few tens of Mq 
similar to that of low luminosity Seyfert galaxies, but ^2 orders of magnitude less than 
that of luminous quasars. Further frequent photometric/spectroscopic monitoring is 
needed to support or refute the different alternatives discussed here on the nature of 
SDSS J153259.96-003944.1. 

Key words: quasars: general - quasars: individual: SDSS J153259. 96— 003944.1 



1 INTRODUCTION 

The diverse observational characteristics of Active Galactic 
Nuclei (AGN) have been reconciled in an unification scheme 
(Antonucci 1993; Urry & Padovani 1995). The basic idea is 
that all the AGNs have broad and narrow line emitting re- 
gions (BLR and NLR respectively) and an obscuring torus. 
Most of the diversity in the observed properties is caused 
by differences in our viewing angles to the torus axis. It 
has also been found that some of the physical parameters of 
AGNs show statistically significant correlations. For exam- 
ple, tight relationships exist between: (i) the radius of the 
BLR and the continuum luminosity (Corbett et al. 2003) 
and (ii) the black hole mass and the velocity dispersion of 
the host galaxy (Onken et al. 2004 and references therein). 
With the advent of new very large surveys, some objects 
that appear to depart from this standard AGN picture are 
beginning to emerge. Understanding these observations are 
important to get a clearer picture of AGN formation and 
evolution. 

Several AGNs were recently found to have peculiar 
emission line characteristics. In the Sloan Digital Sky Sur- 
vey (SDSS), a few high redshift quasars have been dis- 
covered without emission lines (Anderson et al. 2001; 



Hall et al. 2004). Among them the z = 4.67 quasar 
SDSS J153259.96-003944.1 (hereafter referred to as SDSS 
J1533— 00) first reported by Fan et al. (1999) is the object 
of interest in the present study. The optical spectrum of this 
quasar is featureless redward of the Lya forest region; blue- 
ward of ~6800 A, the spectrum has features due to Lya 
absorption at z — 4.52. Near-IR observations confirm the 
presence of a point source and an extended nebulosity in 
SDSS J1533-00 (see Hutchings 2003). Another object (2QZ 
J215454. 3— 305654) with similar characteristics is also re- 
ported by Londish et al. (2004). 

The weak (or vanishing) emission line spectrum seen in 
SDSS J1533— 00 is typical of BL Lac objects, where the lines 
are presumably swamped by the beamed and boosted con- 
tinuum (Urry & Padovani 1995). BL Lac objects are char- 
acterized by strong radio and X-ray emission, optical vari- 
ability and strong (and variable) optical polarization due 
to synchrotron radiation from a relativistic jet. Apart from 
the featureless spectrum, SDSS J1533— 00 does not possess 
these other characteristic properties of BL Lac objects. It 
was not detected in deep radio observations (3<r upper limit 
of 60 ^iJy) and was not found to be optically polarized (with 
a 3(J upper limit of 4%; Fan et al. 1999). It was found to be 
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Table 1. Log of observations and the results of photometry. Here, 
QSO— A, QSO— B and A— B refer to the differential instrumental 
magnitudes respectively between the quasar and the comparison 
stars A and B, and between the comparison stars themselves 



Date 



A Exp. 
Time 
(sec) 



QSO-A 

(mag) 



QSO-B 

(mag) 



A-B 

(mag) 



17.01.00 


R 


1200 


1.74 


± 


0.06 


2.93 


± 


0.05 


1.19 


± 


0.02 




I 


1800 


0.94 


± 


0.05 


2.01 


± 


0.05 


1.07 


± 


0.02 


28.03.00 


R 


1800 


1.92 


± 


0.06 


3.12 


± 


0.06 


1.20 


± 


0.02 




I 


1800 


1.27 


± 


0.05 


2.31 


± 


0.04 


1.04 


± 


0.03 


30.03.00 


R 


1800 


1.79 


± 


0.06 


2.99 


± 


0.05 


1.20 


± 


0.01 




I 


1800 


1.24 


± 


0.05 


2.29 


± 


0.05 


1.05 


± 


0.02 


09.04.00 


R 


1800 


1.81 


± 


0.05 


3.00 


± 


0.05 


1.19 


± 


0.01 




I 


1800 


1.25 


± 


0.06 


2.34 


± 


0.05 


1.09 


± 


0.02 


01.05.00 


R 


1800 


1.82 


± 


0.05 


2.98 


± 


0.05 


1.17 


± 


0.02 




I 


1800 


1.24 


± 


0.08 


2.32 


± 


0.07 


1.08 


± 


0.02 


21.01.01 


R 


1200 


1.97 


± 


0.06 


3.13 


± 


0.06 


1.16 


± 


0.02 




I 


1200 


1.36 


± 


0.07 


2.35 


± 


0.07 


0.99 


± 


0.02 


26.01.01 


R 


900 


2.10 


± 


0.09 


3.27 


± 


0.09 


1.18 


± 


0.02 




I 


900 


1.24 


± 


0.08 


2.25 


± 


0.08 


1.01 


± 


0.04 


19.02.01 


R 


1800 


1.99 


± 


0.04 


3.20 


± 


0.05 


1.21 


± 


0.03 




I 


1800 


1.43 


± 


0.04 


2.46 


± 


0.05 


1.03 


± 


0.03 


25.03.01 


R 


300 


2.05 


± 


0.11 


3.21 


± 


0.11 


1.16 


± 


0.02 




I 


300 


1.32 


± 


0.07 


2.33 


± 


0.07 


1.02 


± 


0.03 


23.04.01 


R 


1200 


2.08 


± 


0.05 


3.25 


± 


0.05 


1.17 


± 


0.02 




I 


1200 


1.42 


± 


0.05 


2.44 


± 


0.05 


1.02 


± 


0.02 



an extremely weak X-ray source and was not detected in 
pointed Chandra observations (Vignali et al. 2001). 

The absence of strong emission lines in the case of SDSS 
J1533— 00 can be due to one of the following reasons: (i) the 
optical continuum being Doppler boosted as in the case of 
BL Lacs; (ii) the optical continuum being amplified due to a 
gravitational microlensing event by a star in an intervening 
galaxy; and (iii) absence of line emitting gas in the vicinity of 
the central UV continuum source. As frequent photometry 
of the source can confirm or reject the first two possibilities 
we have carried out photometric monitoring observations. In 
section 2 we outline the observational program, data reduc- 
tion procedure and the results of the photometric monitor- 
ing. Section 3 discusses the nature of the source. Finally, our 
conclusions are given in Section 4. The cosmological model 
we consider in our study is Q m = 0.3, Qa = 0.7 and Ho = 
70 km s" 1 Mpc" 1 . 



2 OPTICAL OBSERVATIONS AND ANALYSIS 
2.1 Photometry 

Optical photometric observations in Cousins R and I band 
were carried out on 10 nights between January 2000 and 
April 2001 using the 104 cm Sampurnanand telescope of 
the Aryabhatta Research Institute of Observational Sciences 
(ARIES), Nainital. This is an RC system with a f/13 beam 
(Sagar 1999). The detector used was a cryogenically cooled 
2048 x 2048 CCD mounted at the cassegrain focus. Each 
pixel of the CCD corresponds to 0.37 (arcsec) 2 and the entire 
CCD covers 13' x 13' on the sky. Observations were done in 
2x2 binned mode to improve S/N. Typical seeing was ~2" 
during most of our observations. 
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Figure 1. Differential light curves (DLCs) of the quasar SDSS 
J1533-00 with respect to two stars A (RA(2000) = 15:33:07.310; 
Dec(2000) = -00:39:04.0) and B (RA(2000) = 15:33:04.110; 
Dec(2000) = -00:40:49.00). The top most panel gives the stan- 
dard R magnitudes of stars A and B. The other panels give the 
DLCs of the pair of comparison stars and between quasar and 
the comparison stars as indicated within the panels. The left and 
right panels refer to R and I bands respectively. The solid lines 
are the linear least squares fits to the data. 



Initial processing of the images (bias subtraction, flat 
fielding and cosmic ray removal) were done using IRAF 1 
routines, whereas PSF fitting photometry was carried out 
using the routines in MIDAS 2 . To look for variability, differ- 
ential light curves (DLCs) of the quasar were generated with 
respect to stars A and B situated on the observed frames. 
Typical error in our photometry is around 0.03 mag for the 
reference stars and between 0.03 and 0.1 mag for the quasar. 
The log of observations and the results of the photometry 
are given in Table 1. The DLCs in R and I filters with re- 
spect to stars A (RA(2000) = 15:33:07.310; Dec(2000) = 
-00:39:04.0) and B (RA(2000) = 15:33:04.110; Dec(2000) = 
—00:40:49.00) are shown in Fig. 1. From Fig. 1, it appears 
that the quasar shows a gradual fading during the period of 
our observations in both R and I bands, though of course 
more complex behavior cannot be excluded. In particular, as 
we have gaps in our lightcurve, any non-linear fluctuations 



The Image Reduction and Analysis Facility (IRAF) is dis- 
tributed by the National Astronomy Observatories, which is oper- 
ated by the Association of Universities for Research in Astronomy 
Inc. (AURA) under cooperative agreement with the National Sci- 
ence Foundation 

2 Munich Image Data Analysis Systems; trademark of the Euro- 
pean Southern Observatory (ESO) 
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Figure 2. Differential R band light curve of the quasar SDSS 
J1533-00 with respect to star A (RA(2000) = 15:33:07.310; 
Dec(2000) = -00:39:04.0) including the observations of Fan et 
al. (1999) (shown as open circles). Continuous and dotted lines 
give respectively the linear least squares fit and lcr error to the 
data. 



(such as flares) cannot be ruled out. Linear least squares fits 
to the QSO DLCs in both R and I bands give a similar slopes 
of ~0.20 mag/year. However, the QSO DLCs in I band show 
large scatter compared to the R band data. This is due to 
the worse quality of the I band data due to problems in 
flat fielding. Also, since the R band is near the maximum 
of the response curve of the CCD, we consider only the R 
band data in any further analysis. T he amplitude of variabil- 
ity was calculated using A max = \J (D max — D min ) 2 — 2a 2 , 
where D max (Dmin) are the maximum (minimum) in the 
quasar DLC and a the average error in the corresponding 
DLC (Romero et al. 1999). The object was found to show a 
peak to peak amplitude of variability (A max ) of about 0.34 
mag between January 2000 and April 2001. 

SDSS J1533— 00 was also observed for two epochs (27 
June 1998 and 21 March 1999) by the SDSS team (Fan et 
al. 1999). No clear variation is observed between these two 
epochs. Converting their i* magnitudes using equations (1) 
and (2) given by Fukugita et al. (1996) and the zero point 
given by Bessel (1979), the R band magnitudes are esti- 
mated as 19.66 and 19.68 on 27 June 1998 and 21 March 
1999 respectively. In order to compare these quasar mag- 
nitudes to its magnitude during the epoch of our observa- 
tions, we converted the quasar's instrumental magnitudes 
during our observations to apparent R magnitudes by the 
technique of differential photometry using Star A at the po- 
sition (RA(2000) = 15:33:07.310; Dec(2000) = -00:39:04.0 
and whose R band magnitude of 18.46 mag is taken from 
USNO-B; Monet et al. 2003) available on the quasar CCD 
frame. The apparent R magnitude of the quasar thus cal- 
culated during our period of observations range from 20.20 
to 20.56 mag. Thus, combining our observations with those 
reported by the SDSS team, we find that the quasar has 



faded by ~0.9 mag over a period of about 3 years from June 
1998 to April 2001. Using a simple least squares fit to all the 
available data we get (see Fig. 2), 

Am(maj) = (1.17 ± 0.07) + (0.34 ± 0.03) x t(yr). (1) 

Here, Am is the differential magnitude of the quasar with 
respect to the reference star A and t is the time in the ob- 
server's frame. 

2.2 Spectroscopy 

We use the optical spectrum of SDSS J1533— 00, obtained 
by Fan et al. (1999) with the Keck II telescope using the Low 
Resolution Imaging Spectrograph, to obtain the optical con- 
tinuum spectral index, a u (defined through, S v oc v a °), and 
an upper limit of the Lya emission line flux. The observed 
spectrum is well fitted with the existence of Lya absorption 
system at z abs = 4.58 with N(Hi) = 1.1 xlO 21 cm -2 , Lya 
emission line with a FWHM of 105 A (or a velocity width 
of 4630 km/s) at an emission redshift of z em = 4.67 and 
ao = —0.8. N(Hl) in the absorbing gas cloud is obtained by 
self-consistently fitting Lya, Ly/3 and Ly7 lines. The fit to 
the spectrum is shown in Fig. [3] Inclusion of Lya emission 
in the fitting improves the fitting of the Lya absorption line. 
However, as noticed by Fan et al. (1999) the higher Lyman 
series lines do not show complete absorption suggesting the 
gas producing the absorption cannot be a strong damped 
Lya system. Thus our fits will over predict the N(Hl) in 
the absorbing gas and the derived upper limit of the Lya 
emission line flux. 

The estimated flux of the fitted Lya emission line is 
^ 1.06 x 10 -15 erg cm~ 2 s~ x . This corresponds to an emitted 
luminosity of ^ 2.2 x 10 44 erg s _1 . Using the fitted contin- 
uum and emission line flux we get an upper limit (~ 9.9 A) 
on the rest frame Lya emission line equivalent width. 

2. 2. 1 Excess continuum flux 

As discussed before, the lack of emission lines in the spec- 
trum of SDSS J1533— 00 can be due to the continuum emis- 
sion being enhanced either by Doppler boosting of the rel- 
ativistic optical jet emission or by gravitational lensing of 
the continuum but not the line emitting BLR. From the 
observed spectrum we obtain an estimate of the enhance- 
ment needed in the optical continuum so that the emission 
lines appear weak. This is done by calculating the equiva- 
lent width of Lya emission line from the observed spectrum 
of SDSS J1533— 00 and comparing it with the typical Lya 
emission line equivalent width observed in the composite 
spectrum of quasars (Francis et al. 1991). The equivalent 
width (EW) is defined as 

W ob3 = —= (2) 

fJLrcont 

where fi is the magnification. The rest frame W b s is 9.9 A. 
When compared with the typical EW of about 80 A(Francis 
et al. 1991) observed in quasars, this gives a magnification 
of ^ 8. Thus roughly an order of magnitude amplification in 
the continuum flux is needed to explain the observed weak- 
ness of the Lya emission line if it is not due to a lack of 
emitting gas. 
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Figure 3. The observed optical spectrum of the quasar using the LRIS spectrograph on Keck II (solid line) along with the fitted 
absorption and emission lines. The dotted line shows the modeled spectrum containing only the Lya emission line whereas the dashed 
line shows the modeled spectrum containing both Lya emission and absorption lines of Lya, Ly/3 and Lyy. 



3 DISCUSSION 

3.1 Nature of optical variability of SDSS 
J1533-00 

The absence of strong emission lines in the spectrum of SDSS 
J1533— 00 is similar to that of BL Lacs. It is also believed 
that the lack of emission lines in BL Lacs is due to the con- 
tinuum being Doppler boosted. One of the defining char- 
acteristics of BL Lacs is that they show large amplitude 
flux variability over the entire electro-magnetic spectrum. 
Here, we discuss the observed variability properties of SDSS 
J1533— 00 and the known properties of the different kinds 
of AGNs derived based on long term optical monitoring to 
understand the nature of SDSS J1533— 00. The R band ob- 
servations correspond to ~1230 A in the QSO rest-frame. 
A linear least squares fitting to our observations (see Eq. 1) 
combined with two other epochs of observations of the SDSS 
team, gives a rate of change of ~0.35 mag/yr in R band in 
the observed frame. This when transformed to the rest frame 
of the quasar (taking into account the effects of time dila- 
tion; trest = t 6 S (l+ z)) gives a rate of decline in the R band 
magnitude of ~1.9 mag/yr. This rate of decline is an aver- 
age, in reality the decline could be more dramatic if there 
had been a flare in the intervening period. 

Comparing SDSS and POSS measurements on a large 
sample of SDSS quasars, de Vries et al. (2003) have reported 
that the long term quasar variability is consistent with a de- 
caying intrinsic light-curve of the form, 1.1 x exp(— t(yr)/2). 
Similar variability time-scales (~1 year) and amplitudes 



(~0.16 mag) have also been found by Trevese et al. (1994) 
from a 15 year monitoring of 35 quasars between 0.6 < z < 
3.1. Giveon et al. (1999) have reported an average rate of 
change in the intrinsic B band magnitude of 0.28 mag/year 
using 7 years of monitoring data on 42 low- luminosity PG 
quasars sample with z < 0.4. For the same data, Cid Fernan- 
des et al. (2000), using structure function analysis, report 
a variability amplitude and time-scale of 0.18 mag and 1.8 
years respectively. Recently, using large database of variabil- 
ity on SDSS quasars, Ivezic et al. (2004) report a variability 
amplitude and time-scale of 0.32 mag and 1 year, respec- 
tively. Thus, it appears that the normal population of QSOs 
seem to show an intrinsic rate of change of < 0.5 mag/yr. 
In contrast, from the published lightcurves of blazars moni- 
tored over a 18 year time baseline (Webb et al. 1988), con- 
sidering the long term trends, we notice that blazars on an 
average show a larger rate of variability (> 1 mag/yr). 

From the above discussions, it is clear that the observed 
R band variability of 1.9 mag/yr shown by SDSS J1533— 00 
is larger than that seen in typical QSO population. As SDSS 
J1533— 00 is highly luminous (M 4r = —26.6 mag), X-ray 
and radio quiet and at high redshift, based on the existing 
correlations (de Vries et al. 2004; Vanden Berk et al. 2004; 
Cid Fernandes et al. 1996; Cristiani et al. 1996; Giveon et 
al. 1999), we expect it to show a lower rate of change of 
magnitude. However, the observed high rate of variability 
is consistent with that seen in BL Lac objects. Thus our 
observations are consistent with SDSS J1533— 00 being a 
low luminosity AGN (with intrinsically weak emission lines) 
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with its continuum being boosted by relativistic beaming. 
However, the continuum emission mechanism in the jet of 
this object needs to be very different from typical BL Lac 
objects, since its broad-band properties are so discordant. 

Although the sparse long-term variability nature of 
SDSS J1533— 00 is similar to BL Lacs, intra-night mon- 
itoring observations too are needed to confirm the pres- 
ence of relativistic optical continuum emitting jet in SDSS 
J1533— 00. This is due to the fact that the observed large 
amplitude and frequent variability of BL Lacs on intra-night 
timescales (Stalin et al. 2004; Gopal-Krishna et al. 2003) is 
now widely believed to be due to the inhomogeneities in the 
outflowing relativistic jet. Such a study is planned in the 
coming observing season. 



3.2 Is variability due to microlensing? 

In this section, we explore the possibility that the observed 
decline in the light curve of SDSS J1533— 00 is due to mi- 
crolensing of the continuum source by an intervening object. 
Such a model was proposed by Ostriker & Vietri (1985) to 
explain BL Lac objects (see also Gopal-Krishna and Sub- 
ramanian 1991, for a relativistically moving source). Mi- 
crolensing in the intervening galaxy is believed to be the 
cause of abnormal emission line ratios seen among differ- 
ent components in the multiply imaged systems such as 
Q2237+030 (Huchra et al. 1985), APM 08279+5255 (Lewis 
et al. 2002) and SDSS J1004+4112 (Richards et al. 2004). As 
the QSO is not multiply imaged, the line of sight to the QSO 
most probably samples the outer region of the intervening 
galaxy. However, we note that the required magnification 
(i.e., 10) of the continuum for the emission lines to be in- 
visible in SDSS J1533— 00 is much higher than that required 
in the case of multiply imaged systems discussed above (i.e., 
~2). 

For the continuum source (and not the broad line emit- 
ting region) to be magnified significantly due to microlens- 
ing, the microlensing length scale, i.e., the Einstein radius, 
should be larger than the projected radius of the continuum 
emitting region (^ 10 14 cm) and smaller than the radius of 
the BLR (few 10 17 cm). The Einstein radius (Re) is defined 
as 

_ jAGM D ls 

where Di, Di s and D s are the lens, source-lens and source 
angular diameter distances, respectively, and M is the mass 
of the lens. For a lens mid-way between the source and the 
observer (which corresponds to a lens redshift of 0.678), Re 
is ~ 0.01 pc for 1 Mq lens. Such a lens will satisfy the re- 
quirements for microlensing described above. We notice that 
to get the required high magnification the impact parameter 
has to be less than 1.5 xlO 15 cm at the lens plane. The rate 
of change of magnification will depend upon the velocity of 
the star for a given impact parameter. 

We obtain the expected light curve in the observer's 
frame assuming the velocity of the lens to be 200 km s _1 in 
the lens plane for different impact parameters (using equa- 
tion 20 of Refsdal 1964). The results are shown in Fig. 0] 
where we plot the log of magnification as a function of 
time. As we do not know the intrinsic flux and the epoch 
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Figure 4. Simulated microlensing lightcurves (log of magnifica- 
tion versus time) for three impact parameters. Our observations 
(filled circles) as well as those of Fan et al. (1999) (open circles) 
are overplotted on the modeled lensing lightcurves. The impact 
parameter in units of 10 15 cm is given against each curve. The 
dotted horizontal line indicates a magnification of 10. 

at which SDSS J1533— 00 was at its maximum magnifica- 
tion, we shifted the observed points to match the predicted 
curve with similar slope. From Fig. 4 it is clear that the 
sparsely sampled observed light curve can match the pre- 
dicted lensing curve for an impact parameter of ~ 7 x 10 14 
cm for an assumed solar mass lens moving at 200 km s . 
The required impact parameter will be slightly higher if we 
assume the velocity of the lensing star to be higher than 200 
km s . This simple exercise demonstrates that microlens- 
ing is a viable option in this case. Interestingly, the QSO is 
found to be surrounded by excess density of galaxies, with 
the closest companion galaxy at a separation of only ^3.5" 
from the QSO (Hutchings 2003). We need closely sampled 
multi-band light curves to confirm or refute this microlens- 
ing scenario. Also, if microlensing is the correct option, then 
we now expect the Lya emission line to become more clearly 
visible in the spectrum of the QSO as the enhanced contin- 
uum continues to decline. Further photometric and spectro- 
scopic observations should thus constrain the microlensing 
hypothesis. 

3.3 How much gas is there in the BLR? 

The upper limit we derive on the Lya line flux (see Sect. 
2.2) can be used to obtain a limit on the mass of the line 
emitting gas in the BLR. If BLR gas emits very efficiently 
then the mass of the BLR (Mblr) is given by 

M BLR = 5.1(10 u /n e )(L L ya/10 45 ) M (4) 

For the upper limit of Lya luminosity of SDSS 1533—00, we 
get Mblr, < 1.1 M for n e = 10 11 cm -3 . As pointed out by 
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Table 2. Results of photoionization model computations. Here, N(Hl) = total hydrogen column density, nn = hydrogen density, ~Rblr 
= radius of the BLR, log U = logarithm of ionization parameter, L(Lyct) = luminosity of L(Lyo) emission line, N c = number of Lya 
emitting clouds and Mblr = mass of the BLR 
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Baldwin et al. (2003) the mass estimated using the above 
equation will be the minimum mass in the BLR. We obtain 
a realistic upper limit to Mblr using the photoionization 
code CLOUDY (Ferland et al. 1998). 

The ionization state of the line emitting gas is quanti- 
fied using a dimensionless ionization parameter, U. This is 
defined as, 

£912 



U 



4:TvR^ LR hvcn 



(5) 



where n, L912 and Rblr are, respectively, the number den- 
sity of the gas particles, the Lyman continuum luminosity, 
and the radius of the BLR around the central engine. We 
obtain L912 = 2.7 X 10 46 erg s _1 using the observed flux 
above the Lya emission line and a spectral index, a a — —0.8 
(see Sect. 2.2). For a standard AGN, Rblr can be estimated 
from the rest-frame luminosity at 5100 A(i.e., L5100) using 

0.68 



Rblr = 27 A 



lO^ergs" 1 



light days, 



(6) 



(see Corbett et al. 2003 and references therein). The above 
relationship is established using reverberation mapping 
studies. Assuming such a relationship holds good even for 
SDSS J1533— 00, we estimate the radius of the BLR to be 
0.8 pc. Thus we have a constraint on U for an assumed value 
of n. For example, when we consider n = 10 9 cm -3 the con- 
sistent value of log[/ is —0.29. 

We then estimated the Lya emission line luminosity 
from a single cloud (L(Lya) cai ) for a given n, U, total hy- 
drogen column density N(Hl) and ionizing spectrum of the 
QSO using the photoionization code CLOUDY (Ferland et 
al. 1998). This together with the observed upper limit on 
the Lya emission line flux (F(Lya) obs ) is used to get the 
number of BLR clouds, 



47rd 2 F(Lya) c 



(7) 



L(Lya) cai 

Here, for simplicity we assume all the BLR clouds to be 
identical. In addition if we assume the clouds to be spheri- 
cal we can get the mass of the Lya emitting clouds in the 
BLR. The total mass of the BLR contained in N c clouds is 
calculated using 

Mblr = N c x ^Tvr 3 P (8) 
o 

where p and r are the total hydrogen mass density and radius 
of the cloud (r = N(Hi)/2n), respectively. 



The ionizing continuum used in the model calculations 
is a combination of a UV bump (assumed to be a black- 
body with a temperature of 10 K) and an X-ray power law 
(Georgantopoulos et al. 2004) of the form /„ oc ^~ 2 '°. The 
UV and the X-ray continuum slopes were combined using 
an UV to X-ray logarithmic spectral slope of a ox = —2.0. 
This value of a ox is consistent with the upper limit derived 
by Vignali et al. (2003). 

For this assumed incident ionizing continuum shape and 
solar chemical abundances, we carried out grids of photoion- 
ization model computations for various values of N(Hl) and 
n. The chosen range in N(Hi) considers BLR to range from 
optically thin to optically thick. The ranges of other chosen 
parameters are also consistent with those frequently used 
in BLR modeling (see Korista et al. 1997). The results of 
these model calculations are summarized in Table 2. It is 
found that the calculated values of the Lya emission line 
luminosity are weakly dependent on the input values of a ox . 
We also notice that changing the temperature of the UV 
bump from 10 5 K to 10 6 K changes the luminosity by less 
than a factor of 2. Assuming the observed continuum is en- 
hanced by a factor of 10 (due to either Doppler boosting or 
gravitational lensing), we also carried out photoionization 
model computations by de-magnifying the measured con- 
tinuum luminosity by a factor of 10. The results are shown 
within brackets in Table 2. 

Clearly our model calculations are consistent with 
Mblr < 50 M Q in the case of SDSS J1533-00. This is 1 to 
2 orders of magnitude less than that derived for the stan- 
dard high luminosity QSOs (Baldwin et al. 2003). However, 
the estimated upper limit is consistent with BLR masses 
computed for low luminosity Seyfert galaxies (see Peterson 
2004). 



4 CONCLUSION 

We have presented optical photometric monitoring of the 
peculiar quasar SDSS J1533— 00 for a duration spanning 
about 500 days during which the object has varied by about 
0.4 mag. These observations, when coupled with two other 
epochs of observations available in the literature, indicate 
that the quasar gradually faded by 0.9 magnitude during 
the period June 1998 to March 2001. This transforms to a 
variability of ~1.9 mag/yr in the rest frame of the quasar. 
Such a large amplitude of variability is similar to the long 
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term variability nature of BL Lacs. Nevertheless, the lack 
of X-ray and radio emission and optical polarization sug- 
gests that the continuum emission in the jet needs to be 
very different from BL Lacs. Available photometric data on 
the source could not rule out microlensing as the cause of 
the variability as well as the observed lack of emission lines. 
Further monitoring observations are needed to constrain the 
microlensing scenario. 

Photoionization model calculations show the BLR mass 
to be consistent with low luminosity Seyferts (Peterson 
2004), but ~2 orders of magnitude lower compared to those 
expected for high luminosity quasars (Baldwin et al. 2003). 
It is also possible that SDSS J1533— 00 belongs to an un- 
known population of highly luminous AGNs without a BLR. 
This would be in line with the "naked" AGNs discussed by 
Hawkins (2004). Furthermore, if the observed lack of emis- 
sion lines in the discovery spectrum of SDSS J1533— 00 (Fan 
et al. 1999) is due to the continuum being amplified either 
due to Doppler boosting or to gravitational lensing, then we 
might expect to see emission lines emerge as our observa- 
tions show the quasar to be in the declining phase. Further 
photometric and spectroscopic observations could clarify the 
peculiar nature of SDSS J1533— 00. 

There now exist quite a few other objects with feature- 
less optical spectra resembling BL Lacs but which lack the 
significant radio and X-ray emission found in BL Lacs (An- 
derson et al. 2001; Hall et al. 2004; Londish et al. 2002). 
Whether SDSS J1533— 00, along with such other objects, 
fit into the orientation based unification scheme, or whether 
they instead belong to a hitherto unrecognized population of 
radio-quiet BL Lacs or lineless radio-quiet quasars remains 
an open question. 
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